Exact wave functions of electron in a quantum dot 
with account of the Rashba spin-orbit interaction 
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We present the exact wave functions and energy levels of electron in a two-dimensional circular 
quantum dot in the presence of the Rashba spin-orbit interaction. The confinement is described by 
the realistic potential well of finite depth. 
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The Rashba spin-orbit interaction in semiconductor 
quantum dots has been the object of many investigations 
in recent years (see [H, [3, |3, |j] and references therein). 
The Rashlaa interaction is of the form 0, Q 



with standard Pauli spin-matrices 
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The Rashba interaction can be strong in semiconductor 
heterostructures and its strength can be controlled by an 
external electric field. 

The quantum dots in semiconductors can be described 
as effectively two-dimensional systems in confining po- 
tential Vc{x, y). A confining potential is usually assumed 
to be symmetric, Vc{x,y) = Vc{p),p = Then 
the single-electron wave functions satisfy the Schrodinger 
equation 



(3) 



where n is the effective electron mass. There are two 
model potentials which are widely employed in this area. 
The first is a harmonic oscillator potential. Such a model 
[1, |3| admits the approximate (not exact) solutions of 
Eq. (3). The second model P, 01 is a circular quantum 
dot with hard walls Vc{p) = for p < po, Vc{p) = oo for 
p > Pq. This model is exactly solvable. In the framework 
of above models the number of allowed energy levels is 
infinite for the fixed total angular momentum. 

In this paper, we examine a model which corresponds 
to a circular quantum dot with a rectangular potential 
well of finite depth: Vdp) = for p < po, Vdp) = Vq > 
for p > Pq . We shall present the exact coordinate 
wave functions and energy levels for the wells of arbi- 
trary depth. 

The analogous model with replacement Vc{p) — > 
Vc{p) — Vo was proposed in jj, where the approximate 



investigation was developed in the momentum represen- 
tation for a shallow well. As to the coordinate repre- 
sentation, unfortunately the paper Q does not contain 
the radial Schrodinger equations as well as the explicit 
and complete expressions for solutions of these equations. 
Emphasizing different aspects of the same problem the 
present paper and the paper [3] mutually complement 
each other. 

The Schrodinger equation is considered in the cylin- 
drical coordinates x = p cos Lp,y — p sin Lp. Further it is 
convenient to employ dimensionless quantities 



2M 2p 2p 2t. 
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As it was shown in [111, Eq. (3) permits the separation of 
variables: 



«'™(r. If) = u{r)e' 



w{r)e 







(5) 



m = 0,±1,±2, . . . 
due to conservation of the total angular momentum Lz + 

n . 

In [1|, the requirements u(l) = w(l) = were im- 
posed. In the present model, we look for the radial wave 
functions u(r) and w{r) regular at the origin r = and 
decreasing at infinity r — > oo. 

We consider two regions r < 1 (region 1) and r > 1 
(region 2) separately. 

In the region 1 (w = 0) the radial equations may be 
written in the suitable form 
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+ r—^ {k^k^r"^ - (m + l)^)wi 



where 
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Note that in the case /3 = {k^ = k^) Eqs. (6) are the 
Bessel equations. Therefore, following [l|, [3] we use the 
known properties [§| 

d n 



J„(fcr) 



-fcj„+i(fcr), 



dr 
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of the Bessel functions in order to obtain the exact solu- 
tions 

ui{m,e, f3,r) = ci/i(m, e, /?, r) 

+ digi{m,e,P,r), 

Wi{m,e,f3,r) = cigi{m + l,e, f3,r) 

+ difi{m + l,e,f3,r) (9) 

of system (6), where 

1 



/i(m,e,/3,r) 



(Jm(fci r) + Jm{kfr)) 



gi{m,e,(3,r) = - (Jm(fci r) - J™(fcj^r)) 



(10) 



are the real linear combinations of the Bessel functions 
with real arguments. Here Ci and di are arbitrary coef- 
ficients. The radial wave functions ui{r) and wi{r) have 
the desirable behavior at the origin. 

In the region 2 (v > 0) we have the radial equations 
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{k+k^ + (m + lf)w2 
= - k^)r^ - ^^2) , 



where 



h (e,w,/3) = W-y 



(11) 
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Now in the case /? = {k^ = /c^) Eqs. (11) are the 
modified Bessel equations. Therefore, using the known 
properties [1| 
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of the modified Bessel functions it is easily to get the 
exact solutions 

U2im,e,v,P,r) = C2/2(to, e, /5, r) 

+ d2g2{m,e,v,(3,r), 

W2{m,e,v, P,r) = C2g2{m + 1, e,v, (3,r) 

- d2f2{m + l,e,v,f3,r) (14) 



of system (11), where 



f2im,e,v,(3,r) = - {K^ik^ r) + K,n{k^r)) , 
52(m,e,t;,/3,r) = {Kraik2r) ^ Kmik+r)) (15) 



are the real linear combinations of the modified Bessel 
functions with complex arguments. Here C2 and d2 are 
arbitrary coefficients. At large values of r the functions 
/2(to, e, D, /3, r) and (72("^, e, D, /3, r) behave as 



/2(m,e,w,/3,r) 



52(m,e,w,/3,r) 
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where 7 is defined as follows. 



- e - ^2/4 . 

cos(7) = -^^ -= , sm(7) = 

\/v — e 
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Thus, the radial wave functions U2{r) and W2{r) have 
the appropriate behavior at infinity. It should be noted 
that the radial wave functions have the infinite number 
of zeros for the finite value of energy. 

Since the solutions (9) are valid for e > —fP/A and 
the solutions (14) are valid for e < w — /3^/4 then the 
complete energy range is —0^/4: < e < v — 0^ / A. 

The continuity conditions 



ui{m,e,P, 1) 
u'i{m, e, (3, 1) 
wi(to, e,/3, 1) 
w[{m,e,/3, 1) 



U2{m,e,v,P, 1), 
U2(m,e,v,/3, 1), 
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(18) 



for the radial wave functions and their derivatives at the 
boundary point r = 1 lead to the algebraic equations 




M{m,e,v,P) 



for coefficients ci,C2,di and d2, where 



(19) 
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TABLE I: Energy levels. 

V P e 



TO = 



25 





3.98 


9.94 


19.61 


25 


1 


3.49 


9.81 


19.15 


25 


5 


-4.40 


2.83 


13.40 



25 10 -23.25 -18.31 -9.67 

49 4.41 11.13 22.75 35.91 

49 1.4 3.49 11.03 21.87 35.79 

49 7 -10.40 -3.71 9.55 24.22 

49 14 -47.11 -41.45 -32.19 -19.60 -2.25 

100 4.77 12.09 24.97 40.08 60.28 81.84 

100 2 2.97 11.75 23.30 39.73 58.65 81.42 

100 10 -21.91 -16.95 -4.88 14.82 35.04 56.69 

100 20 -97.96 -91.83 -81.75 -67.58 -49.91 -28.53 -1.66 

m = 1 

25 9.94 17.46 

25 1 9.02 17.85 

25 5 -2.52 9.26 

25 10 -23.22 -17.29 -4.85 

49 11.13 19.85 35.91 

49 1.4 9.41 20.50 34.28 

49 7 -9.62 1.84 20.95 36.73 

49 14 -47.04 -41.12 -31.51 -13.33 

100 12.09 21.66 40.08 57.25 81.84 

100 2 8.85 22.59 37.08 58.13 79.02 

100 10 -22.91 -14.79 4.24 31.61 55.83 74.93 

100 20 -97.91 -91.72 -81.26 -66.83 -48.24 -17.86 

TO- = 2 

25 17.46 

25 1 16.13 

25 5 1.30 16.08 

25 10 -22.86 -13.35 -0.20 

49 19.85 30.35 

49 1.4 17.37 31.72 

49 7 -6.88 10.11 32.20 

49 14 -44.83 -40.77 -26.05 -4.15 

100 21.66 33.34 57.25 76.20 

100 2 17.08 35.51 52.95 78.21 

100 10 -22.12 -8.70 16.99 49.86 74.91 

100 21) -l)7.«l -<)i.:i9 -81). 2!) -(io..-)2 -.-iT.O!) 15 



M{m,e,v.p) = 



( /i(m,e,/3,l) -f2{m,e,v,P,l) gi(m,e,/3,l) -g2{m,e,v, (3,1)\ 

f[{m,,e,0,l) -f!^{m,e,v,/3,l) g[{m,, e, (3,1) -.g^(m, e, t;, /3, 1) 

gi(m + l,e,/3, 1) -g2{m + l,e,v,P,l) /i(to + 1, e, /3, 1) f2{m + l,e,v, (3,1) 

\g[{m + l,e,p,l) -ff^(m + 1, e, t;, 1) /{(m + 1, e, /3, 1) /^(m + 1, e,t;,/3, 1) / 

I 



(20) 



Hence, the exact equation for energy spectrum is It should be stressed that in the explored model the num- 

det M{m,e,v, P) =0. (21) 
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FIG. 1: Radial wave functions for m = 1. Solid line for u{r), 
dashed line for w{r). 
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FIG. 2: Radial wave functions for m = —2. Solid line for 
u(r-), dashed line for w(r). 

ber of admissible energy levels is finite for the fixed total 
angular momentum. Note that Eq. (21) is invariant un- 
der two replacements m — > —(to + 1) and /? —[3. 
If the energy values are found from Eq. (21) then it 



is simply to get the values of coefficients ci, C2, di and d2 
from Eq. (19) and the following normalization condition 
{v?{r) + w^(r)) rdr = 1. 

Now we present some numerical and graphic illustra- 
tions in addition to the analytic results. 

Table U shows the dependence of energy e on the 
Rashba parameter /3, the well depth v and angular mo- 
mentum number to. The values of (3 are 0, 0.2-^1;, ^Jv and 
2^/v for the given value of v. We see that the number of 
energy levels decreases if the well depth v decreases and 
if angular momentum number to increases. Besides we 
see that the negative energy values appear for the posi- 
tive V when the values of [3 becomes sufficiently large in 
comparison with the value of ^/v. If /3^/4 > v then the 
energy spectrum is completely shifted to negative region. 

Figures [1] and [5] demonstrate the continuous radial 
wave functions for to = 1 and to = — 2 respectively 
while /3 = 2,1; = 100 and e = 37.0825 in both cases. 
The solid lines correspond to the functions u{r) and 
the dashed lines correspond to the functions w{r). The 
numerical values of coefficients are ci = 4.22035, C2 = 
-4067.87, di = -0.7139284 and = 880.843 in the case 
of Fig. 1 and ci = 0.713928, C2 = 880.843, di = -4.22035 
and d2 = 4067.87 in the case of Fig. 2. We see that the 
radial wave functions rapidly decrease outside the well. 
Figures indicate some symmetry properties of the radial 
wave functions under replacement to — > —(to -1-1). 

In our opinion the examined exactly solvable model 
with the realistic potential well of finite depth is physi- 
cally adequate in order to describe the behavior of elec- 
tron in a semiconductor quantum dot with account of the 
Rashba spin-orbit interaction. Further we intend to gen- 
eralize our consideration by including the magnetic field 
effects on the orbital motion of electron in a quantum 
dot. 

The author thanks K. Pankrashkin for pointing out the 
reference f7l|. 
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